Vortex mechanics in planar nano-magnets 
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Collective- variable approach for study of non-linear dynamics of magnetic textures in planar nano- 
magnets is proposed. The variables are just arbitrary parameters (complex or real) in the specified 
analytical function of complex variable, describing the texture in motion. Starting with such a 
function, a formal procedure is outlined, allowing to obtain the (non-linear) system of differential 
equations of motion for the variables. The resulting equations are equivalent to Landau-Lifshitz- 
Gilbert dynamics as much as definition of collective variables allows it. Apart from collective- variable 
specification, the procedure does not involve any additional assumptions (such as translational in- 
variance or steady-state motion). As an example, equations of weakly non-linear motion of magnetic 
vortex are derived and solved analytically. A simple formula for dependence of vortex precession fre- 
quency on its amplitude is derived. The results are verified against special cases from the literature 
and agree quantitatively with experiments and simulations. 
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Unlike bulk magnets or thin films (whose magnetic tex- 
ture consists of domains and domain walls, extensively 
studied during the last century in the framework of mi- 
cromagnetics) nano-magnets, due to their restricted ge- 
ometry, display markedly different textures, consisting of 
vortices and anti- vortices, some of which can be bound to 
the surface This opens up a wealth of possibilities 

for new physics, not studied in the established science of 
micromagnetics. Some micromagnetic approaches (such 
as Thiele equation^ for steady-state uniform transla- 
tion of magnetic texture) can be successfully applied to 
study of motion of single magnetic vortex in bigger nano- 
magnets 0,0], but the full account for restricted geometry 
(their most essential feature) and multi- vortex dynamics 
requires a different formalism, which is free from assump- 
tions, specific to micromagnetics. 

Nevertheless, Thiele equation is still the most use- 
ful analytical tool for study of nano-magnet dynamics. 
Another approach, based on volume averaging of the 
Landau-Lifshitz-Gilbert (LLG) equation in vector form, 
produces qualitative agreement with micromagnetic sim- 
ulations and the experiment 0, H ■ But it leads to un- 
derestimate of texture mobility [8| and vortex precession 
frequency 0, because LLG equation is non- linear (due to 
constraint on the length of magnetization vector, masked 
when it is written in vector form) and it's volume aver- 
aging (a form of linear superposition) is, generally, not 
justified. Consideration of spin- waves on magnetic vortex 
background [ij. [loj does reproduce translation of magnetic 
vortex and predicts higher-energy spin-wave modes, but 
is limited to linear consideration of small vortex displace- 
ments only. Including higher order terms in deviation of 
magnetization from (magnetic vortex) background is not 
only mathematically hard, but also bound to have dif- 
ficulties reproducing the complicated non-linear motion 
of the multi-vortex texture, which may completely de- 
part from the original static background. Non-uniform 



background also makes it difficult to deal with non-local 
dipolar forces, which is the reason why in many of such 
works (Refs[^ andfiol in particular) dipolar interaction is 
replaced by local in-plane anisotropy, making the equa- 
tions partial differential (instead of integral partial dif- 
ferential). This is justified in the limit of vanishing nano- 
element's thickness, but quantitative agreement with ex- 
periments and simulations in wide range of geometries 
is only possible when magnetostatic interaction is fully 
accounted for. 

Here another approach to magnetization dynamics in 
planar nano-dots is proposed. It is a collective-variable 
theory, capable of dealing with complex multi- vortex con- 
figurations (fully describing relative motion and deforma- 
tion of constituent vortices). Its only approximation lies 
in the definition of collective variables, which are just 
an arbitrary parameters in a complex function of com- 
plex variable. Given such a function, the approach pro- 
duces a system of ordinary differential equations (ODEs) 
of motion (with no integral terms even in the presence 
of magnetostatic interaction) for these variables. They 
are exactly equivalent to Landau-Lifshitz-Gilbert equa- 
tion, apart from the restriction, implied by the collec- 
tive variable definition. The approach assumes neither 
translational invariance, nor steady state motion and is 
fully capable of describing non-linear dynamics (if one 
can solve the resulting non- linear ODEs). The external 
field, spin-torque and other potential energy terms can 
be easily added without sacrificing simplicity (dynamical 
equations still remain ODEs). It also allows to include 
phenomenological dissipation, akin to Gilbert's damping 
term in LLG. As an illustration, the equations of motion 
for linear and weakly non-linear magnetic vortex dynam- 
ics in circular nano-dots are derived and solved. As a 
check, the well known analytical result for vortex preces- 
sion frequency (originally obtained by solving the Thiele 
equation) is then recovered in the limiting case of large 
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dots. 

The original motivation for this work comes from re- 
cently published description of low-energy (single- and 
multi-vortex) magnetic configurations in planar nano- 
dots of arbitrary shape in terms of functions of complex 
variable The present approach can be thought as a 
way to "animate" these configurations by making them 
move in accordance with LLG dynamics. Despite this, 
one may easily generalize it to other sets of trial func- 
tions without substantial sacrifice in simplicity. 

The usual starting point for consideration of magne- 
tization dynamics is LLG equation It is very well 
suited for numerical computations, but is not good for 
analytical ones. This is because the non-uniform effec- 
tive field in LLG, around which magnetization vectors 
precess, depends, in turn, on the whole magnetization 
vector field (if dipolar interaction is properly taken into 
account). This makes it non-linear integral partial differ- 
ential equation, which is extremely hard to solve analyt- 
ically. Therefore, instead of LLG, let us go deeper and 
consider, as a starting point, kinetic Lagrangian density, 
which was first introduced by Doring 12] : 
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where 9 and ip are polar and azimuthal angle of the 
magnetic moment M in a spherical coordinate system, 
t is time, 7 ~ 1.76 • 10 11 rad/(sT) is gyromagnetic ra- 
tio, Mg is saturation magnetization and 9q is a constant. 
Parametrization of M via spherical angles conveniently 
satisfies the constraint \M\ = Ms, leaving only two of 
its components independent (and bounded). The sys- 
tem of two Euler-Lagrange equations for the extremum 
of the corresponding action over 9 and tp with additional 
time-independent potential energy terms subtracted, are 
equivalent to Landau-Lifshitz eauation[l2l [l3| . They do 
not depend on 9q, which can be used to ensure that r is 
zero at the boundary of the magnet. 

Collective variable approach is then similar to the Ritz 
method [3] of solving boundary value problems: first, 
one selects a trial function, parametrizing a wide set of 
possible solutions, and then finds the values of parame- 
ters, giving the correct answer (extremalizing a certain 
functional, as per variational principle). Ritz method in 
its original formulation finds wide applications in micro- 
magnetics for solving static problems. Dynamics is not 
much different. One may look for the extremum of action 
of the full Lagrangian, including kinetic and potential 
energy, parametrized by a certain set of scalar parame- 
ters. Condition for this extremum produces dynamical 
equations for the parameters, allowing to compute their 
evolution in time. 

While the above general recipe is applicable to arbi- 
trary choice of trial functions, to make further consid- 
eration more specific, let us focus on a particular very 



general family Consider a cylinder, made of soft fer- 
romagnetic material, placed in the origin of the Carte- 
sian coordinate system in such a way that the axis Z 
is perpendicular to the cylinder face X>, which is not 
necessarily circular, distribution of magnetization vec- 
tor M(r), r = {X,Y,Z} inside can be assumed to be 
independent on the coordinate Z . It can be conveniently 
parametrized by a complex function w(z,z) of complex 
variable z = X+iY (line on top denotes the complex con- 
jugation, so that z — X — iY), and to express normalized 
magnetization m = M /Ms as 

2^ 

m x + im Y = — — (2) 
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m z = ±— — =, (3) 
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which automatically satisfies the constraint m — 1 . The 
sign in mz controls polarization of the vortex core: mz — 
±1 at the vortex center. 

In flat cylinder the equilibrium static magnetization 
distributions can be represented as a combination of soli- 
ton and meronQ 

r m/ci \f(z)\ < ci 

H*>*) = \ f{z)/\jf{z)~m ci < \f(z)\ < c 2 , (4) 
(f(z)/c 2 |/(*)|>ca, 

where f(z) is an analytic (df(z)/d~z = 0) function of 
complex variable z and c\, ci are real constants. 

Besides the complex variable z the function w(z, z) 
usually depends on a number of other variables like c±, 
c-2, and others hidden inside f(z). The simplest case of a 
single magnetic vortex (Usov's anzats[lj) corresponds to 

i(z-A) 
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with c\ = 1 (absorbed into vortex core radius ry) and 
C2 — > 00 in (j4]) , since there are no anti- vortices. This mag- 
netization distribution depends on real parameter ry and 
complex parameter A, which are the collective variables 
in this case. The parameter A — ax +iay is the position 
of the vortex center. The problem, considered below, is 
how to find dynamical equations for these (and other sim- 
ilar) collective variables, assuming they are ry = ry(t), 
A = A(t) functions of time t. 

The Landau-Lifshitz equation can also be written di- 
rectly in complex notation [15j|. However, as discussed 
earlier, our starting point will be the kinetic part of the 
magnetic Lagrangian density ([1]), expressed through the 
collective variables and their time derivatives. In the 
complex notation its ingredients are 
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so that (p} can be rewritten as 
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where dot over the variable denotes the time derivative 
and w(z, z, t) from depends on time t via the collective 
variables in the trial function f(z). The meron part of 
dU), as per our selection of 9q = 7r /2, gives no contribution 
to the kinetic Lagrangian, because \w(z, ~z)\ — 1 inside it, 
and therefore r ~ mz ~ cos = 0. 

To derive dynamical equations total Lagrangian is 
needed, which is the Lagrangian density, integrated over 
the particle volume 
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where Lz is cylinder's thickness and T>s is part of cylin- 
der's face, occupied by soliton @, for which \f(z)\ < c\ 
or 1/(^)1 > c-i- It can be simplified for arbitrary / in (U) 
by noting that 
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1 - w(t)w(t') 1 dw(t) 
1 + w(t)w(t') w(t) dt 



(11) 

where the variable t' is considered independent and does 
not take part in differentiation. Interchanging the Im op- 
eration and time derivative with the area integral (which 
is possible because the area element is real and integra- 
tion is a linear operation) we arrive at the following ex- 
pression for the total kinetic Lagrangian 
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where for the purpose of calculating the area integral 
and differentiating it is assumed that all the collective 
variables inside w and definition of integration region T>' s 
(note the prime) depend on the new independent time 
variable t', whereas inside w they still depend on t (only 
after differentiation the t' is replaced by t). It was also 
noted that inside the soliton T>s the function w is analytic 
and does not depend on z, while its conjugate w does not 
depend on z. This formula can be checked directly. It 
can be further simplified if there are no boundary-bound 
vortices and anti- vortices. In this case it is possible to 
integrate (I12[) by parts, making use of Greene's formula 
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for any reasonably good function u, which yields 
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where function arguments are omitted, but it is still as- 
sumed that all the collective variables inside w depend 
on t, but all of them inside w and definition of the region 
T>' s depend on t' . 

The expression (TT^|) (and (TT4"|) for the case of solitons 
fully contained inside the particle) is the main result of 
this paper. To derive equations of motion for collective 
variables Xi(t), entering the trial function /({a;,}, z), it is 
now sufficient to write the full Lagrangian: 

L({i,}, {x t }) = T({xi}, {x t }) - U({xi}), (15) 

where U is the potential energy (including exchange, 
magnetostatic and, possibly, other energy terms); and 
use it to derive the system of Euler-Lagrange equations 
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(16) 



which extremize the corresponding action. This allows 
to treat problems of (multi-) vortex dynamics as prob- 
lems of classical mechanics. It is also worth noting that 
apart from restrictions, implied by a particular choice of 
/, defining collective variables Xi, the above considera- 
tion involves no approximations and corresponds to the 
solution of LLG equation exactly. 

Let us now proceed with examples. 

First, to illustrate that, similarly to classical mechan- 
ics, this theory may suffer from poor selection of trial 
functions, consider dynamics of a uniformly displaced 
magnetic vortex ([5]) with ry = const and A = A(t). One 
may readily check that all three expressions for kinetic 
Lagrangian (fTU]) . (fT^]) and (fT4")l yield in this case. It 
means that the Euler-Lagrange equations (|16p reduce to 
conditions of static equilibrium, or (in the case of poten- 
tial energy, consisting of the exchange and magnetostatic 
terms) that such undeformed vortex stays in the center 
of the cylinder. This is similar to the conclusion from 
micromagnetics that moving walls always have their pro- 
file different from that of stationary ones[l||. In some 
sense, one may say that modification of magnetic tex- 
ture (of domain wall or a magnetic vortex) is the way 
how it "remembers" that it is moving. 

Let us now turn to a more complex trial function, de- 
scribing the vortex displacement without formation of 
magnetic charges on cylinder's side [2], which is particu- 
lar case of a more general class of trial functions]^] 
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where again ry = const and A — A(t) = ax(t) + iay(t), 
\A\ < 1/4 is a complex pair of collective coordinates. The 
vortex center, where f{zc) = is located not exactly at 
Zc = A, as in the case of uniformly displaced vortex, 
but rather at z c = (1 - \] \ - AAA/ (2 A). While it is 
possible to change variables and write the equations of 
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motion for zq (t) directly, let us illustrate one of the pow- 
ers of the present approach, which is a great freedom in 
selecting the parametrization, and write them for A(t). 
Also note that here the parameter a in the original ex- 
pression of Ref. is substituted by — lA, which changes 
nothing but makes the vortex center displacements coin- 
cide in phase with the complex parameter A. That is, 
real A correspond now to real Zc- This again is a matter 
of convenience and does not change anything, since the 
parametrization can be arbitrary. The total Lagrangian 
up to the second order in |A| from (|T4"f and (|T5|) with full 
account for vortex core shape deformation is 

C = ±k 2 {a x {t)a' Y {t)~a Y {t)a' x {t)) - k 2 (a 2 x (t)+a 2 Y (t)) , 

(18) 

where C = L/([iQMgirLzR 2 ), &2 are dimensionless and 
k 2 = (1 + r v (4 log2 — 3))/ (jfi Ms) has units of seconds; 
k 2 is second order expansion coefficient of the poten- 
tial energy (which includes exchange and dipolar inter- 
action). The equations of motion (fTBf arc 



n 2 a' x (t) ± k 2 ay{t) = 
K2a!y{t) =F k 2 a x (t) = 



(19) 



which, for initial conditions a x (0) — ao, ^y(v) — 0, have 
the following solution 



ax{t) = ao cos(wf) 
ay{t) = ±aosin(wi), 



(20) 



corresponding to the circular motion of the vortex around 
the dot center with frequency uj = ujq = k 2 /n 2 . It is im- 
portant to note that the direction of vortex motion is 
not arbitrary. It depends on its core polarization, but 
not chirality, since T and U are independent on the sign 
of w or /. The vortices with mz = — 1 at the center 
rotate clockwise, while the vortices with mz — 1 coun- 
terclockwise. This is in full agreement with simulations 
and experiments of Ref. 1171 but in disagreement with its 
conclusions, since vortex chirality (included in "handed- 
ness" ) plays no role in determining the direction of vor- 
tex rotation. This also allows to guess that the vortex 
core polarization in the simulation of Ref. was posi- 
tive, which is natural to assume, but was not specified 
by the authors. A similar polarization of the core can 
be guessed from Fig. 2 in Ref. 0, but with significant 
uncertainty, since it is masked by low resolution of the 
measurement in Y direction, as discussed therein. 

To make a more rigorous quantitative confirmation of 
the present theory, let us compute rotation frequency of 
magnetic vortex in the limit of large flat circular dots 
with Lz <C R and R ^> Le, where Le = ^C/fioM'g is 
the exchange length, C is the exchange stiffness, and R is 
dot's radius. In this case rv <C 1, and k 2 ~ l/(jfioMs). 
The second order expansion of the energy of the vor- 
tex ([TT]) with vortex core neglected, was first published 
in Ref. [Ill- If the exchange contribution of the order 



Le /R -C 1 is neglected in that expansion, the coefficient 
k 2 in large dots is fully determined by the energy of vol- 
ume magnetic charges |l8[. Converting to SI units, for 
precession frequency we get 
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where /ms(i) = 1 — (1 — e~ x )/x, g = Lz/R and 
G ~ 0.915966 is Catalan's constant. This expression 
(apart from measurement units and the value of numer- 
ical constant in (1221) . which is exact here) coincides with 
the expression for vortex frequency, obtained in Ref. @ on 
the basis of Thiele equation and quantitatively confirmed 
there by experiments on large dots. 

Real magnets also inevitably dissipate energy of mov- 
ing spins in the form of heat. But Lagrangian formalism 
in its pure form does not include dissipation. It is added 
externally via a Rayleigh dissipation function D 



D 



(23) 



which is then included into the right hand side of the 
Euler-Lagrange equations (1161) via an additional term 
—dD/d£i. The matrix Dij consists of phenomenolog- 
ical dissipation coefficients. Judgin g fr om the abstract 
of the unpublished report by Gilbert [19( it is possible to 
speculate that Lagrangian formalism was also his start- 
ing point and his dissipative term (whose full microscopic 
justification is still an open problem [20j]) has similar ori- 
gins. Thus, must be related to Gilbert's phenomeno- 
logical dissipation constant. This relation is, probably, 
best established by considering the energy balance in the 
system, but let us leave it for now as an open problem 
and treat as independent phenomenological param- 
eters. A choice of D 
solution (l20l) into 



d(a' x (t) 2 + a' Y {t) 2 ) changes the 



a x {t) = ao e- dk2t/{d2+ ^ ) cos (k 2 tn 2 /{d 2 + n 2 2 )) 
a Y (t) = ±a e- dk2t/(d2+K ^ sin (k 2 tn 2 /{d 2 + k 2 2 )) .(24) 

As one can see, similarly to the case of linear oscilla- 
tor, the vortex precession frequency starts to depend 
(slightly) on (small) damping coefficient. 

Finally, let us consider weakly non-linear vortex dy- 
namics by taking into account kinetic and potential en- 
ergy terms, corresponding to the 4-th order in \A\. Con- 
tinuing the expansion of kinetic Lagrangian (|14[) with the 
trial function (TIT)) leads to the following expression 

C = ±(k 2 + Ki{a 2 x {t) + a 2 Y {t))) (a x (t)a Y (t)-a Y (t)a' x (t)) 
-k 2 (a 2 x {t) + a 2 Y {t)) - h (a 2 x (t)+a 2 Y (t)f , (25) 
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where k 4 = 2 - 4(23 + r v { (6061- 6397r^)/8- 1152(1- 
r v ) log 2))/(jfj, Ms) like k 2 has units of seconds and &4 is 
the next potential energy expansion coefficient. The cor- 
responding equations of motion become non-linear, but 
they are solved exactly by (f2"0")l with 
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plitude via potential energy expansion coefficients. Ex- 
pressions for kinetic Lagrangian in (TT8|) . (|25|) for the trial 
function (|17[) can be reused in other calculations, includ- 
ing time-dependent external field, spin-torque and other 
potential energy terms. One may expect them to be as 
simple as the examples above. 

I'd like to thank Vladimir N. Krivoruchko for reading 
the manuscript and many valuable suggestions. 



Similarly to other non-linear oscillators, the vortex rota- 
tion frequency becomes dependent on the rotation am- 
plitude. Derivation of the expressions for hz and k 4 in 
general case with full account for vortex core deforma- 
tion is rather cumbersome and, together with the anal- 
ysis of their dependence on the dot dimensions, will be 
the subject of another forthcoming paper. Nevertheless, 
preliminary versions of these expressions are attached in 
the form of Mathematica file as a supplemental material. 
They can be used to compute vortex precession frequen- 
cies for various dot geometries, not covered here. 

Concluding, several equivalent alternative expressions 
for kinetic Lagrangian ©, (ITUl) . (fT2"]) . (|14|) of an arbitrary 
trial function /, defining collective variables in (possibly 
multi- vortex) magnetic texture Q in flat nano-elements 
are derived. They allow to obtain non-linear equations of 
motion for these variables similarly to the ones in classi- 
cal Lagrangian mechanics. Apart from collective variable 
definition, this theory is exact and involves no additional 
approximations beyond those of Landau-Lifshitz-Gilbert 
equation. It is validated here by considering magnetic 
vortex precession in cylindrical nano-dot. In the limit of 
large flat dots its frequency coincides with experimental 
data and known theoretical estimations, based on Thiele 
equation^. The question of the direction of vortex ro- 
tation is elucidated. It is found to depend on vortex 
core polarization only and not on its chirality. Also, 
the analytical solutions for vortex rotation in dissipative 
magnet are derived, its frequency (|24|) is dependent on 
phenomenological damping parameter. Finally, weakly 
non-linear rotation of vortex is considered, allowing to 
establish the relation (|26|) between its frequency and am- 
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